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Abstract 

In this paper, we present the theoretical studies of a refractive index map to implement a Gauss 
to Jo-Bessel-Gauss convertor. We theoretically demonstrate the viability of a device that could be 
fabricated on Si/ SiiyOy/ SiixyGcxCy platform or by photorefractive media. The proposed device 
is 200//m in length and 25/im in width, and its refractive index varies in controllable steps across 
the light propagation and transversal directions. The computed conversion efficiency and loss are 
90%, and —0A57dB, respectively. The theoretical results, obtained from the beam conversion 
efficiency, self-regeneration, and propagation through an opaque obstruction; demonstrate that 
a 2D graded index map of the refractive index can be used to transform a Gauss beam into a 
Jo-Bessel-Gauss beam. To the best of our knowledge, this is the first demonstration of such beam 
transformation by means of a 2D index-mapping which is fully integrable in silicon photonics based 
planar lightwave circuits (PLC). The concept device is significant for the eventual development 
of a new array of technologies, such as micro optical tweezers, optical traps, beam reshaping and 
non-linear beam diode lasers. 
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I. INTRODUCTION 



In 1987 Durnin introduced non-diffracting Jp Bessel beams as members of a family of 
solutions to the homogeneous Helmholtz equation, with distinctive properties, e.g. their 
transverse profile does not change as the beam propagates in free space and exhibit self- 
regeneration when disturbed by non-transparent obstacles [Ij. Among them the zero order, 
Jo, Bessel beam has captured most attention due to its peculiar intensity distribution, which 
is focused on the propagation axis, besides non-diffracting propagation 

Whilst mathematically viable, Bessel beams are not square integrable, i.e. the beam 
contains infinite energy, thus, rendering the ideal solution infeasible ^|4j. Yet good approx- 
imation exists, namely Bessel-Gauss beams. Introduced by Gorin et al. [4] these solutions to 
Helmholtz homogeneous equation resemble the ideal Bessel beam with the sole difference 
that they bear finite power; and, therefore, can be realized experimentally, while retain- 
ing the primordial characteristics of self-reconstruction and diffraction-free propagation for 
lengths of interest to many optical applications, e.g. optical tweezers and microscopy [IHTj. 

Experimentally, Durnin et al. [8] were the first to show that it is possible to generate 
Bessel-Gauss beams by means of a circular slit and a lens placed one focal length away. 
Hakola et al. [9] proposed a Nd:YAG cavity with a diffractive mirror to transform the Gauss 
beam from a pump laser into a Bessel-Gauss beam of arbitrary order; Arrizon et al. [10] and 
Otero [llj independently used holograms to alter Gauss beams and produce zero and first 
order Bessel-Gauss beams; lately, Zhan [12j has shown that it is possible to generate them 
by using a radially polarized beam and surface plasmon resonance. Alternatively Cong [2] 
used phase elements to generate zero order Bessel-Gauss beams. Recently novel designs have 
used waveguides to produced diffraction-free beams, v.gr. Canning [13j presented a Fresnel 
waveguide as a diffraction-free mode generator, and Ilchenko et al. [14J experimentally 
showed that cylindrical waveguides can produce truncated Bessel-Gauss beams; and Tsai 
et al. showed that is possible to use an acoustic tunable lens to generate this family of 
beams. 

The afore mentioned methods produces high quality beams. Nonetheless their design 
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and fabrication on photonic integrated circuits (PIC) or planar lightwave circuits (PLC) 
would involve complex processing stpdf, rendering the integration into most PICs platforms 
impractical. It is a matter of debate whether resonator/ waveguide designs (see for example 
[T3l [T4] ) could be used for this purpose. In our experience on such devices, the quality of 
the beam, and conversion efficiency critically depend on the coupling efficiency. Moreover 
the converted beam free space propagation is limited to ^ 1mm [14j. Hence the motivation 
to attain a fully integrated Gauss to Jo-Bessel-Gauss micro-convertor based on current Si 
foundry technology is justified. If viable, it could be readily integrated into semiconductor- 
chip-based optical elements of great importance to a manifold of optical applications. For 
instance they have been shown to increase the scanning resolution and tissue penetration 
depth of optical coherence tomography (OCT) systems by ^ 50% [16j . 



II. PROBLEM DESCRIPTION AND MODELING 

In Durnin's primordial design, an impinging Gauss beam diffracts from the surface of an 
axicon lens P, [U [TT] , forming a converging wavefront and thus giving raise to the Bessel- 
Gauss beam. Analytically it can be described by means of diffraction at the Fresnel limit. 
However accurate this method could be, it presumes the outbound field to extent infinitely 
in the transversal direction in order to comply with the boundary conditions [6j. To avoid 
this complication, we will assume hereafter the output beam to be Jo-Bessel-Gauss as 
defined by Gori et al [4J, i.e. the superposition of a Gauss and Jo-Bessel profiles. 

It could be contended that only a PIC-integrable-axicon lens would suffice to achieve 
the desired conversion. Such lens, though, is equivalent to the fabrication of a continuous 
alteration of the index of refraction, albeit homogeneous, over a conic volume, posing practi- 
cal and non-trivial integration problems when it is manufactured using the well-established 
silicon fabrication technologies, namely the difficulty to shape the conic surface. To alleviate 
this conundrum, we seek a practical proffle of the refractive index that can act as feasible 
substitute for the axicon lens, suited for integration into current silicon foundry and manu- 
facturing technologies. 
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Before proceeding any further, it is important to remark that non-singular graded refrac- 
tive indexes (GRIN) cannot yield perfect non-diffracting beams as proven by Hayata fT8] . 
However a refractive index mapping can potentially achieve high quality conversion from 
a Gauss beam to a Jo-Bessel-Gauss beam, which is different from the fundamental Bessel 
beam presented by Durnin et al ^y. 

In Figure [l} the input and output beams are known a priori^ and hence the problem is 
reduced to determining the heterogeneous refractive index in region A, which will result 
in the transformation of the input to output fields in region B. This is the optical inverse 
scattering problem [T9H22]. 



Analytically, the input and output fields are described by harmonic functions of the 
form Ei = 0(A;,r)e^^^ and Eo = C{k^f)e'^^'^^ respectively; where (j) and ( are solutions to 
the homogeneous Helmholtz equation. In the transformation region, medium A, the field 
is Ea = 'd{k^r)r{t)] where d is determined by the optical properties of the medium. If the 
material is linear, then r(t) = e^^^, and the problem narrows to determining the shape of 
the permittivity, 6:, and permeability, /x, that realize the transformation (f) ^ after some 
finite propagation distance, d. 



In reality, the manipulation of permeability, /i, is complex. On that account the transfor- 
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Figure 1: Schematic representation of a Gauss to Jo-Bessel-Gauss transformation via a 
heterogeneous medium, the transformation device in ragion A has a refractive index map 
that needs to be calculated in order to echieve the desired beam transformation. The 
device is embedded in a homogeneous medium 5, where Bessel and Gauss beams are 
known solutions to the Helmholtz equation. 
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mation medium A is set to be non-magnetic, with some constant value for //. Consequently, 
Maxwell's equation for the electric field E for this region is written as: 

V^^-A^£(r)^ = -V0^.V£), (1) 

where 6:(r), the permittivity, is a real value function and r = (u^v^w) describes a position 
vector in the transformation space. Outside the device the material is homogeneous, i.e. 
£{f) = constant for all r G S, where both Gaussian and Bessel functions are well known 
solutions to equation [ij Inside the device, the permittivity is limited to the range achievable 
in practical fabrication methods: 

emin < e{f) < Smax, M T ^ A. (2) 

Recalling that if Ea = '^{r)r{t)^ then equation [l] can be reduced to the time independent 
Helmholtz equation: 




Observe that if the field ??(r) is given, then the only unknown in equation [s] would be 
6:(r), which is the permittivity map that we wish to determine. Henceforth let us examine 
^(f). 

Denote by Si and S'o, the boundary surfaces between region A, and the incident Bi and 
exit Bo surfaces, respectively. Then ?9(r) can be split into the three regions: 

10(f) V feSi, 
i^{f) V feA, (4) 
C(f) V feSo, 

where : A ^ is the transforming field. Since we wish to control and thence decrease 
the complexity of the problem, we define the shape of across region A. Accordingly, 
we introduce the real valued functions f^g:A^^ [0, 1] and j : A ^ i?^, continuous 
and square integrable, i.e. smooth functions. These functions are invoked to describe the 
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transforming beam as a superposition of the desired input and output fields, in addition to 
an arbitrary field 7 that can be thought to represent the residual beam. The field ip{f) can 
be written now as: 



= + 9{r)C{f) + 7(f), 



(5) 



to comply with the boundary condition f{Si) = g{So) = 1, f{So) = g{Si) = and 7(5^^) = 0. 
Then substituting equation [5] into equation [3] results in a second order partial differential 
equation: 



[f{f^m + gi^a^) -^v • [fim^ + gi^a^) <^ = v [[f{f)m + gi^a^) ■ vin^ 

(6) 

Recalling vector calculus identities, the right hand side of equation |6] can be reduced to: 



f{r)(t){r) + g{r)C{r)) • Vln^(f) = (7/; • V) Vln5(f) + (Vln^(f) •?/;). (7) 



Observe that in equation |6| w and ij) are vector functions with known values across the 
three domains previously described (see equation ^ . Thus by means of equation [t] we can 
rewrite equation |6] as a diflFerential equation that describes the permittivity, 5, as: 

uj^li^e{f) + v) Vln5(f) + (vin5(f) -A^w. (8) 



This is the nonlinear equation for the permittivity that we ought to solve within the 
boundary conditions and fabrication limitations defined earlier: constant permeability, /x, 
refractive index range 1.5 < n(r) < 3.5, where n(r) = ^^(r)/x; and fabrication minimum 
feature size d. 



The computation of an inverse problem as the one hitherto described by equation [8] is not 
trivial, and analytical solutions are rarely found. Moreover, the problem described by the 
above equation is non-linear and ill-defined, i.e. its solution is not unique and is dependent 
continuously on the data [231 [24] . As a result, many numerical schemes have been devel- 
oped to elucidate the solution, v.gr. the adjoint problem [HI [24], linear sampling [251427] . 
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variational methods [28] , and lately transformation optics [22l EHJ [30] . To find a solution 
to this quandary, we developed a numerical technique based on the variational method and 
the adjoint problem. 

We begin by noting that if the input and output beams are symmetrical around the 
propagation axis and does as well the field in the transformation region, then we can scale 
down the problem to two-dimensions: z and the longitudinal and transversal direction, 
respectively. Concerning our matter at hand both Bessel and Jo-Bessel-Gauss beams are 
symmetric around the propagation axis; moreover the transformation field as described by 
the above equations is also symmetric around the axis z. Therefore, our problem can be 
scaled down to a two dimensional numerical scheme. 

Accordingly, the simulation space is discretized on a grid size by dP . The solver is 
given an initial guess of the permittivity map Sij across the region A at every point Zi^yy 
The solution to equation [3] is computed via a variational approach as described in [28] : at 
every point Zi the transversal function of is determined by solving equation [8] to the 
boundary conditions (see equation [5]) : if{y^Zii) and and the field at that point 

ilj{y^Zi). This iterative approach is carried for every point Zi. Thus the transversal shape 

f[r) 




Transversal direction x(a.u.) Normalized propagation direction (r) (a,u.) 

(a) Transformation transition (b) Characteristic functions / and g 

Figure 2: Illustration of (a) superposition of transition stages in the transformation of 
Gauss to Jo-Bessel-Gauss beams; and (b) the transformation functions / and g. 
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of £ij^ and consequently that of n^j, is determined to match the resulting field to the ideal 
transformation given by equation [5| 

This method solely does not guarantee that the solution will be exempted of extreme 
or singular values, withal the variational method is especially sensitive numerical approxi- 
mations [28]. To account for the latter and work out the extrema in the refractive index, 
we implement a solution coming from the adjoint problem [24] . The extreme and singular 
values are approximated to the closest value in the defined range. We then propagate the 
initial Gauss beam through the device and determine how much the field at every point 
diverges from the ideal transformation of equation [5| including as well the near and far 
fields. An iterative refinement on this grid regions follows, as described in [IHlEl], to match 
the numerical computation output to the ideal Jo-Bessel-Gauss beam. The beam propaga- 
tion method is used to propagate the initial Gauss beam through the heterogeneous region A. 

It is important to note that both, the variational and adjoint problem methods, while 
widely tested may render inexact solutions due to numerical approximations which can lead 
to instability and are in general related to the limited resolution of any numerical method. 
To counterweight this difficulties the general approach is to use smaller grid sizes even at 
the expense of computational time increasing exponentially. 



III. RESULTS. 



The transformation functions / and g are plotted in Figure [2b| The initial guess for the 
permittivity map is a Gaussian profile in the transversal direction y for all points along 
width a = 5 and a maximum peak of e^nax = 9 or n^nax = 3. The initial beam in region Bi 
follows a Gaussian function, with width lO/xm and normalized energy to unity. 

The outcome of the calculation and optimization process described earlier is depicted 
in the schematized Figure |3] and the fine resolution (d = Inm) plot of the refractive index 
in Figure [4} The refractive index map is symmetric around the propagation axis z, and is 
comprised of a smooth variation of the refractive index between 1 and ^ 2.57. The minimum 



feature size, d = = (P^ is Inm, and has the final dimensions 200/xm length by 20/xm wide. 
Notice that the minimum feature size defines the smoothness of the refractive index map. 



The initial impinging Gauss beam in region Bi propagates through the device in region 



and it is transformed into a Jo-Bessel-Gauss beam with a 95% efficiency. Figure 5a 



shows the results for the field at 5/xm from the exit of the device, with the corresponding 



far field depicted in Figure [5bj As described earlier, this is the optimized value considering 
the current state-of-the-art nano-fabrication techniques and facilities available to realize the 
device. 



In Figure [5| it is observed that both near and far fields follow a Jo-Bessel-Gauss profile 
with 98% accuracy. The converted beam divergence is less than 1% after propagating 50/im 
and less than 5% after traveling 1mm. The energy density at the principal peak is 95% 
of the converted beam's total energy, which makes it ideal for OCT applications and free 
space light-based communication. Losses due to material absorption, scattering and back 



y 

z 




200 am 



Figure 3: Sketch of the proposed refractive index map for the convent or as defined by the 
transformation functions / and g] with device length of 200/xm and width 25/xm. Albeit 

here shown in 11 steps, the refractive index varies smoothly between 1 to 2.57 in 
controllable steps. The device is rotationally symmetric around the propagation axis z. 
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reflections are 10% of the incident fleld, or 0A57dB^ which can be signiflcantly reduced by 
using antireflection coatings at the input and output surfaces. 



Within the known hmitations in fabrication processes, it is noted that not aU values 
of the index of refraction can be achieved. A theoretical study is required to define the 
tolerance to minimum feature size and its repercussions in the beam profile. We study this 
case by further modifying the grid size, i.e. the minimum feature size, d = = (P from 
Inm to lOnm. The results are shown in Figure [6| 



As it is evident from Figure [6a[ modifying the refractive index grid size results in a beam 
that also resembles a Jo-Bessel-Gauss beam with a 95% accuracy, with the sole difference 
that this beam as it propagates will evolve into a non-Bessel-Gauss beam. This can be 



drawn from the far field, plotted in Figure 6b, Observe that at infinity the field no longer 




200 fim 



Figure 4: Fine resolution refractive index map for the convent or, with minimum feature 
size d = Inm, as defined by the transformation functions / and g. Device length is 200/xm 
and width 25/xm. The refractive index varies smoothly between 1 to 2.57 in 50 controllable 

steps. A practical device could be fabricated on a Si/ SiiyOy/ SiixyGcxCy slab and the 
refractive index steps achieved by controlled oxidation or with photo-refractive materials. 
The device is rotationally symmetric around the propagation axis z. 
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corresponds to the far field of a Jo-Bessel-Gauss beam. Nevertheless, a coarse mapping will 
produce a highly focused beam, where 95% of the energy is localized at the first peak, and 
the beam diffraction over 1mm propagation is less than 5%. 

Bessel-Gauss beams, alike Bessel beams, exhibit self regeneration when their path is 

perturbed by a non- transparent scatterers. To examine this self-healing property of the 

outbound Jo-Bessel-Gauss beam, we insert a random index medium on its propagation 
path. The results are shown in Figure [7| 
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Figure 5: Transverse electric field intensity of the converted beam: (a) simulated near field 
computed at 5/xm from the output of the device (green), and fit, 95% confidence, to 
Jo-Bessel-Gauss (red). Its far field profile (b), and (c) zoom to principal peak. 
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The generated Jo-Bessel-Gauss beam is made to impinge on a random index media, 
which consists of a circle-hke region of 4/xm, set at 50/xm from the output of the device, 
and refractive index of n = 2.0. The scatterer partiaUy obstructs the propagation of the 
beam. The scattered beam profile is computed at 50/xm from the ahen obstruction (see 



Figure 7a and its magnified version in Figure 7c), as weU as the far field (see Figure urn). As 



1.0 - 
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-• -- Simulated near field (modified grid} 
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Figure 6: Transverse electric field intensity profile for the converted beam for two diflFerent 
grid sizes: (a) shows the optimal grid size of Inm (green) and modified grid size of lOnm 
(blue), both computed at the near fields, 5/xm from the output of the device. The far fields 
for both output beams (b), for optimal grid size (green) and modified grid (blue); and (c) 
magnification of the near field around the principal peak. 
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can be seen from the results in Figure [7| the beam described earher (see Figure [5]) exhibits 
self-regeneration. The divergence, computed by comparing the width of the principal peak, 
from the scattered field to the original transformed beam is 0.3%. Since the beam is not a 
full Bessel beam, some light is scattered by the object, as expected, giving rise to a loss of 
19% or OmSdB. 

The 2D mapping for beam conversion heretofore described while challenging to fabricate. 
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(c) 

Figure 7: Transverse electric field intensity profile for: (a) the converted beam (green), its 

fit (red), and the scattered field profile (blue), computed at 50/xm from the scattering 
region. The far field for the undisturbed and perturbed beams are shown in (b); and (c) 
magnification of the fields around the principal peak. 
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could be built using current nano-manufacturing techniques. A prototype, verbi gratia^ 
could be obtain by controlling the degree of oxidation in Si/ SiiyOy/ SiixyGcxCy [32l |33] 
systems, or by horizontally stacking nano-layers of photo-refractive materials such as chalco- 
genide glasses which have been shown to have a wide range of refractive index [34l [35] . The 
advantages of these materials are that they are relatively well-developed, and the material 
systems are compatible with PICs and PLCs. 



IV. CONCLUSIONS 



We have, to the best of our knowledge, theoretically demonstrated the viability of the 
first PIC and PLC compatible device utilizing a heterogeneous refractive index map, to 
achieve the transformation of a Gauss beam into a Jo-Bessel-Gauss beam. The computed 
device has a loss of 0.457rf5, and high energy focus, where 95% of the output beam energy 
is concentrated at the principal peak of the Jo-Bessel-Gauss beam profile. The beam 
has a divergence of < 5% over a propagation distance of 1mm, and exhibits self-healing 
when partially obstructed by an opaque object. The resulting beam has a Jo-Bessel-Gauss 
profile for both near- and far-fields. It is note that our design is based on current manu- 
facturing techniques, such as controlled oxidation of Si/ SiiyOy/ SiixyGcxCy or by stacking 
different layers of chalcogenide glasses. An integrated Gauss to Jo-Bessel-Gauss beam con- 
vertor could enable on-chip applications which take advantage of the beams self-healing and 
non-diffr active properties, such as micro optical tweezers, traps and couplers, photonic inte- 
grated circuits for telecommunications and computerized micro tomography and microscopy. 
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